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Abstract 

The q 2 dependence of the ir-r] mixing amplitude is examined with the use 
of QCD sum rules. The linear slope of the mixing function 9(q 2 ) is found 
to be much smaller than that for p-to mixing. Thus the mixing amplitude is 
approximately the same in the space-like region as in the time-like one, and 
one may neglect the q 2 dependence of the mixing. A comparison between a 
hadron-meson, an effective chiral model, and the QCD sum rules approaches 
is made. 
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I. INTRODUCTION 



The recent interest in examining the off-shell behavior of p-u mixing |I]-|3J] and tt-i] 
mixing is based on the observation of Goldman, Henderson and Thomas (GHT) Jl]] 
that the meson mixing previously used in the calculation of the nucleon-nucleon (N-N) 
charge-symmetry-violating (CSV) potential || may not be correct. In particular, the value 
of p-u mixing found experimentally in the time-like region (q 2 > 0) may not apply in the 
space-like region, from which the CSV potential is generated. Subsequently, other authors 
have discussed this issue for p-u mixing by using different effective models and also 

find appreciable q 2 dependence of the mixing amplitude. With the advent of QCD sum 
rules |J and their phenomenological success in describing hadronic properties [p~Q| p~2|] , it is 
natural to apply this technique to meson-mixing, particularly since the method is closer to 
QCD than other models that have been used. 

In the pioneering work of Shifman, Vainshtein and Zakharov (SVZ) ||, the "on shell" 
p-u mixing was used as an illustration of the power of QCD sum rules. However, in their 
paper some fine points remain elusive and there is no discussion of the "off-shell" behavior 
of meson-mixing. These details were re-examined by Hatsuda, Henley, Meissner and Krein 
(HHMK) [|J, who used both Borel and finite-energy sum rules, together with dispersion 
relations to determine the momentum-dependence of p-u mixing. They found a rapid q 2 
variation of the mixing parameter 6(q 2 ). In this paper we shall follow the spirit of HHMK 
in applying it to tt-t] mixing. By performing the Borel analysis, we can extract 9{q 2 ) as a 
function of the CSV parameters, and we can also study how the QCD condensates affect 
the mixing function 6(q 2 ). 

This paper is organized as follows: In section I, we establish notations and definitions. 
The calculations are described in section III. In section IV, we relate our results to those 
of GHT |ij and Piekarewicz- Williams (PW) 0. Finally, in section V, conclusions and a 
summary are presented. 

II. tt-t] MIXING (FORMALISM) 

The 7r-T] mixing function 8(q 2 ) is defined by the mixed correlator of tt° and 77° 
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vr^(g 2 ) = i J (fx e iqx < T<K°{x)rf(ti) > (1) 

where 7r° and r]° are eigenstates of an isosymmetric strong-interaction Hamiltonian. 

We can express the physical n-r] fields, which are eigenstates of the full strong interaction 
Hamiltonian, in terms of 7T°, f]°, by introducing the mixing matrix 
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where e is a "mixing angle". By saturating the Im tc 7111 with the 7r and t] poles and using a 
dispersion relation for Re ^ v {q 2 ), we can write 

tit?/ 2\ 9(<1 2 ) /o\ 

where the mixing function is related to the mixing angle e as 

6(q 2 ) = [e(mj) - e(m 2 )] g 2 + [mje^ 2 ) - m^e^ 2 )] . (4) 

Since the e "mixing angle" may be a function of q 2 , the value of eim^) may be different 
from that of e(m 2 ); in this way the mixing function 8(q 2 ) develops a linear dependence on 

On the other hand, in order to study the meson properties from the quark-gluon degrees 
of freedom without being plagued by the problem of wave functions, we need to introduce 
interpolating quark currents for the correlation function, and the meson properties are rep- 
resented by the experimentally measurable matrix elements of these quark currents. In the 
present case, we choose to work with the axial vector currents rather than pseudoscalar ones 
for a number of reasons. Foremost is that, although 7r and rj are pseudoscalar particles, the 
use of the axial vector currents gives a better convergence property • Secondly, on the 

mass shell, the two currents are directly related. In an SU(3) notation, we use 

Ju( x ) = -^j^lnlbu{x) + d^ 5 d(x) - 2s7 M 7 5 s(x)) . (5) 



The correlator of these two currents is defined by 

n;2 = i Jd i xe^<Tjl(x)f u (0)> . (6) 



Using the definition of the mixing angle e, (Eq. (3)), the PCAC relation for tc and rj, 

< °l4(°)l 7r (p) >= 1 < Q \j*(°)\v(p) >= * fr, P», where f n and /„ are the decay 

constants for 7r and rj, and the pole approximation for ImnJ^, we have: 
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(7) 



Because of the structure of PCAC and the ir-i] pole approximation, there is no q 2 g iiv 
term in Im ir*%. In addition to the n and rj poles, let us also include the next resonances in 
the Im 7r^; these resonances are a pair of pseudo-vector particles the &i(1235) (denoted by 
A) and /ii(1170) (denoted by B); these resonances contribute to Im ir*%, 



(q»Qv - q 2 qnu)hfB [e'(m 2 A )5(q 2 - m\) - e(m 2 B )5(q 2 - m|) 

where /a, /b are the decay constants for A and B particles. 
Combining these two contributions, we can rewrite 
1 



(8) 



-ImitH 

7T 



q^q u Imn^q 2 ) - q 2 g fiu Imir 2 (q 2 ), 
Imit^q 2 ) = gj(q 2 - ml) - g v 5(q 2 - mj) 

+hA${q 2 - m 2 A ) - h B S(q 2 - m 2 B ) 
+ higher resonances + continuum 
Jm7r 2 (g 2 ) = h>A5(q 2 — m 2 A ) — /ib^(<? 2 — m 2 B ) + higher resonances + continuum , 
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where 



g~K = Uf v e(ml) 

h A = — n— e (m A ) 
mi 



9v = f^f v e ( m l 
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These parameters g n , g v , h A and hs are unknown and will be determined from the QCD 
sum rules. Once they are determined in terms of the CSV parameters and the Wilson 
coefficients of the OPE for the quark current- current correlation function, we can substitute 
these numbers into the mixing function 9(q 2 ) and discuss its momentum dependence in the 
region of interest. 

In the case of ti-t] mixing, in contrast to the case of p-u mixing, there is no "on shell" 
value for 9(q 2 ), due to the large difference of it and rj masses, Am = m v — m T = 412 



e'(m|) 



(10) 



MeV ~ 3m 7T ; this mass difference presumably comes primarily from the s-quark contribu- 
tion in the 7/. There is, thus, no direct experimental method to study tt-t] mixing, unlike 
e + e~ — > 7T + 7r _ at the cj-resonance fll3]]14}| , since 7r and 77 cannot be formed from one photon 
states in e + e~ collisions; however, the decay of the rf to 7r + 7r~7r° via the 7r° pole serves to 
set an "experimental" mixing at the rj mass [[[JJ. Here we shall content ourselves to use 
experimentally measured masses for the ?t(135), 7/(547), b(1235) and /i 1 (1170) particles as 
input parameters, and distinguish between e(m 2 ), e(m 2 ), e'(m^), and e'{m 2 B ) as indepen- 
dent variables. Although we do not rely on an expansion in terms of the mass difference 
5m 2 = m 2 — m 2 to perform the calculations as in the p-uo case, the large mass difference of 
7r and rj is comfortably accommodated without difficulty 



III. CALCULATIONS 

One important ingredient of the QCD sum rules is the operator product expansion (OPE) 
|T6|j for the correlation functions. In our case, we calculate the OPE of ir** by keeping 
operators up to dimension 6. For the Wilson coefficients, we work to first order in a s , a and 
where a s (a) is the fine structure constant for the strong (electromagnetic) interaction, 
and m q is the current mass of quark q (m u for u quark, and for d quark). We use the 
definition Q 2 = —q 2 ; we do not include the q 2 dependence of a s but choose its value at 1 
GeV 2 , a s (1 GeV 2 ) ~ 0.5. Our calculation is similar to that for the p-u case but we are 
dealing with axial vector currents for pseudoscalar mesons; thus, the reader may convince 
himself or herself of the signs of our results by counting the number of 7 matrices in each 
diagram. In the following, we shall neglect a detailed discussion of each diagram (c.f. ref. 
0]) and simply write down the relevant answers: 

^(9) = <?/,<?^i(<? 2 ) - Q 2 9^ 2 (q 2 ), Q 2 = ~q 2 
W^^.nQ^ + l + l+oJi,) 
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where 



d = + 0{m, 



C*2 = 2[m u < uu > —rrid < dd >] 

C 3 = vra s [< uu > 2 — < d<i > 2 ] 

81 
ss 

2 „ 1 j ^ 2 



vra[4 < uu > 2 - < dd > 2 ]. (12) 
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It is helpful to note three points: 

(1) For vector mesons, the conservation of the vector current leads to a ir^ that must be 
proportional to q^q v — q 2 g fiu ] therefore 7i"i = 7r 2 . This is no longer true in the n-7] case, since 
the axial vector current is not conserved. A side benefit is that we now have two independent 
sum rules for 7r^, and we can take advantage of this property to extract information on the 
parameters appearing in the mixing function 8(q 2 ) 

(2) If we restrict ourselves to operators up to dimension 6 and make the vacuum saturation 
assumption (VSA) for the four quark condensate |§, we find that the s quark in the r\ 
current (Eq. (5)) does not contribute at all. The reason is that the s quark could only arise 
in the four quark condensate < qTqsTs > (q = up or down quark) which vanishes under the 
VSA. 

(3) Taking the difference of the sum rules for 7Ti and 7T2 we see that the coefficient C3 for the 
dimension 6, four-quark condensate drops out. Therefore the uncertainty in this coefficient 
due to the assumption of the VSA can be ignored. 

The next step is to relate these CSV parameters and condensates to the phenomenological 
parameters g n , g v , Jia and lis in the RHS of Im tt^ (our model for Im ir^). In order 
to assure convergence, to accentuate the lower resonances and to lessen the effect of the 
continuum, we follow the standard practice of applying a Borel transformation to iri and 7^2 



tarn - ja j^iQ'r . (13) 

fixed 

In this way, we obtain a set of simultaneously linear equations for g n , g v , Ha and Kb'- 
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The sum rule for 7Ti gives 

(m 1 ) + 2" " U/ li 



Co + C 2 (-?-) + l -d (-L 
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+ higher resonances (7r', t]')(A', B') + continuum , 



(14) 



and that for tt 2 yields 



2 vm 6 ; 

4v/3 



M 2 



+ higher resonances(A', 5') + continuum . (15) 



Taking the difference, we get 



2C : 



/J_\ = 4^3 
VMV M 2 



+ higher resonances (ir', rf) . 



(16) 



It is worthwhile to call attention to the fact that the contributions of higher resonances 
6i (1235) and /ii(1170) and continuum have been completely canceled in Eq. (16). The same 
is true for the electromagnetic continuum, whereas higher resonances in the pseudoscalar 
channel, such as 7r'(1300) - ?/(1290), will survive (in principle) in Eq. (16). However, it 
will turn out that taking these resonances explicitly into account does not change the Borel 
analysis noticeably. This is in contrast to the case of p-u mixing where the corresponding 
high resonances (p'-ou') are important. 

By taking the derivative with respect to (jj^z) on both sides of Eq. (16), we obtain 
another sum rule: 

-c 2 



2VS 



= ml g w - m\ g v 



(17) 



With these two equations, we can solve for g n and g n : 



9tt 



Co 



K 2^M 2 ) 

' _c 1 _ 

Finally, the mixing function is 



9r, 



m 2 v /M 2 



' M 2 



mz 



m 2 — m 2 



M 2 



rat 
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6(q 2 ) = a(M 2 )q 2 + (3(M 2 ) , 



(18) 



(19) 



where 



«( M 2 ) EE 9V 9 * 



p 

m , } s mfr.-mfo (2Q) 

where we assume the SU(3) symmetry for the decay constants f n = f v . 

In Fig. 1 we plot the dependence of g n and g v on M 2 . In the preferred region of 1.0 
GeV 2 < M < 1.6 GeV 2 , the difference between g n and g v is small. Furthermore, it can be 
seen from Fig. 1 and Eq. (18) that if M 2 S> m 2 , g n ~ g v . We thus find that 6(q 2 ) has a 
weak q 2 dependence. 

In order to fix the values of g n , g v , a and (3, we choose a suitable Borel window and 
average the functions g 7T (M 2 ) } g^M 2 ) over the range of M 2 (1 GeV 2 ~ 1.5 GeV 2 ). At this 
point, one should note that if we included the higher resonances 7r(1300) and 7/(1295) in 
the phenomenological side of the sum rules and left their masses as parameters, we could 
perform the stability analysis to obtain the "optimal" values for g n and g v . However, our 
numerical results do not suggest such "optimal" values are well-defined. The numerical 
values we obtain for these parameters are 

g n =lU MeV 2 o„ = 117 MeV 2 

y yv (21) 

a = 3.65 x 10" 4 /3 = 3660 MeV 2 

All numbers have an error of about 10-20%, which is due to the uncertainty of the isospin 
symmetry violating quantities md ~ mu and 7 = — 1 in the coefficient C2 and the fact 

771 c i -\- 111, 11 <^WU,J> 

that M 2 > 0.9 GeV 2 for stability. 



IV. COMPARISONS WITH OTHER MODELS 

In this section, we shall compare our results with those of other authors, who use a 
different approach. 

TT-rj mixing has been considered previously by Maltman and Goldman (MG) |5|,[|] and 
Piekarewicz (P) 0. Malt man used chiral perturbation theory || to calculate the n-r] mixing 
to the one- loop order. The mixing function 8(q 2 ) is expressed in terms of meson masses and 
other physical observables; theoretical uncertainty comes from the electromagnetic mass 
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difference of kaons. In order to verify the assumption which underlies the GHT calculation. 
Maltman and Goldman |J also use a chiral quark model to calculate ir-t] mixing and get a 
consistent result with chiral perturbation theory. 

Piekarewicz |7j used a purely hadronic model to calculate tt-tj mixing, where the mixing 
is generated by a NN loop. In order to fix the renormalization point for the divergent 
integral, he chose the "on shell" value of m v to fix the intercept of the mixing function. 

He also derived a ratio between the slope at the origin of the p-u and tt-t] mixing ampli- 
tudes 

< n\H\r] > slope at q 2 = ^ g^g v 

< p\H\u > slope at q 2 = f p g w 

which he claimed to be less mo del- dependent and to hold approximately for a reasonable 
range of q 2 . 

In view of the lack of "on-shell" experimental value for the ir-r] mixing amplitude (the 
previous calculations are based on SU(3) mass splitting |l7| , |T8H of pseudoscalar mesons and 
a pole model analysis of the rj-rj 1 system []T7|,[nj), our calculation provides an independent 
result for the slope and intercept for the mixing function. We list the different results from 
the three calculations in Table 1. 

It is to be noticed that our slopes are three times smaller than that of chiral perturbation 
theory and correspondingly so is the ratio a/ (3. Because of the small slope, our mixing 
function only changes about 5% from q 2 = m 2 to q 2 = —m 2 . Therefore, the mixing function 
is practically a constant. On the other hand, the intercept of our mixing function (/?) is fairly 
close to that obtained in chiral perturbation theory. Furthermore it should be noted that 
both in the GHT and in the PW approaches the q 2 variation of the tt-t] mixing amplitude is 
much smaller than the one for the p-u mixing amplitude obtained in the same approaches 
|l|-|3|]. Finally, our mixing function gives an "on shell" value (q 2 = m 2 ) of 3800 MeV 2 , which 
is slightly smaller (but within our estimated error) than that of the other two approaches 
and the one of 4200 MeV 2 obtained from the CSB NN force (I7HI3- 
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V. SUMMARY AND CONCLUSIONS 



In this paper, we use QCD sum rules to study rr-r] mixing, choosing a combination of the 
two sum rules coming from the correlation functions of two axial vector currents. We are able 
to extract the leading behavior of the mixing function 9 as a function of q 2 without the need 
to use four-quark condensates. We obtain 9(q 2 = m 2 ) = 3800 MeV 2 , and a (^-dependence 
of 0(q 2 ) that is compatible with zero. This is in contrast to the p-oj mixing amplitude, which 
varies strongly with q 2 . Our results are in qualitative agreement with those obtained from 
various quark and hadronic models as well as from chiral perturbation theory. 
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FIGURES 

FIG. 1. The mixing parameters g w and g n as functions of the squared Borel mass M 2 
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TABLES 



TABLE I. Comparisons of different approaches for the tt-tj mixing function 

a (MeV 2 ) a/P (MeV^ 2 ) 6(q 2 = m 2 )(MeV 2 

Chiral Perturbation 1.08 xl0~ 3 3808 2.836 xl0~ 7 4131.68 

Theory |J2| 



4200 



Hadronic Model g 1.7x 10~ 3 3800 4.4x 10~ 7 

QCD Sum Rules 3.5 xlO" 4 3700 9.45 xl0~ 8 3800 
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